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Interference effects strongly affect the transport characteristics of a benzene single-electron tran- 
sistor (SET) and for this reason we call it interference SET (I-SET). We focus on the effects of 
degeneracies between many-body states of the isolated benzene. We show that the particular cur- 
rent blocking and selective conductance suppression occurring in the benzene I-SET are due to 
interference effects between the orbitally degenerate states. Further we study the impact of re- 
duced symmetry due to anchor groups or potential drop over the molecule. We identify in the 
quasi- degeneracy of the involved molecular states the necessary condition for the robustness of the 
results. 

PACS numbers: 85.65.+h, 85.85.+j, 73.63.b 
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I. INTRODUCTION 

Molecular electronics, due to perfect reproducibility 
and versatile chemical tailoring of its basic components, 
represents one of the most promising answers to the in- 
creasing miniaturization demand of information technol- 
ogy. A crucial issue in molecular electronics is thus the 
understanding of the conduction characteristics through 
single molecules 1 -. 

Single-molecule-transport measurements rely on 
the fabrication of a nanogap between source and 
drain electrodes and the formation of a stable 
molecule-electrode contact. Nanogaps are nowa- 
days routinely obtained using different techniques 
including electromigratio n 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 !, mechani- 
cal break-junctio n 11 ! 12 ' 13 ! 14 and scanning tunnelling 
microscop y 15 ' 16 ! 17 . Also the challenging goal of effec- 
tively gating a nanometer sized molecule in presence of 
macroscopic metallic leads has been achieve d 6 ' 14 ' 18 . 

A stable contact between molecule and leads is com- 
monly realized with the mediation of anchor groups at- 
tached to the molecule during its chemical synthesis. 
Also direct coupling of the molecule to the electric leads, 
though, has been very recently reportecU 3 -. One of the 
advantages of the first connecting method is some con- 
trol over the contact configuration of the molecule^ and 
the possibility to design the strength of the tunnelling 
coupling by choosing specific anchor groups 6 - ' 17 ' 20 ' 21 . All 
previous achievements combined with experience accu- 
mulated with semiconducting and carbon-based single 
electron transistors (SET) allowed in recent years to mea- 
sure stability diagrams of single molecule transistor de- 
vices thus realizing molecular spectroscopy via transport 
experiments 2 -' 3 .'!' 5 .'^' 7 .'^'!!!. 

Single molecule transistors display transport proper- 
ties which are very different from those of conventional 
single electron transistors. In fact, vibrational or tor- 
sional modes 7 -! 1 - - and intrinsic symmetries of the molecule 
can hinder or favor transport through the molecular SET, 
visible e.g. in the absence or presence of specific excita- 
tion lines in the stability diagram, or in negative differen- 
tial conductance features. Many-body phenomena as e.g. 
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FIG. 1: (color online) Schematic representation of the two 
different setups for the benzene I-SET considered in this pa- 
per. The molecule, lying on a dielectric substrate, is weakly 
contacted to source and drain leads as well as capacitively 
gated. 



the Kondo effect, have been observed as well 2 -' 3 -' 5 - ' 10 ' 22 . 

Despite the experimental progress, the theoretical un- 
derstanding of the properties of single organic molecules 
coupled to electrodes is far from being satisfactory. 
On the one hand, numerical approaches to transport 
based on the combination of Green's function meth- 
ods with tight binding model or density functional the- 
ory have become standard in the study of transport at 
the nanoscalei. These methods are appropriate to in- 
vestigate quantum transport through molecular bridges 
strongly coupled to leads. In this regime various groups 
have recently discussed the possibility to observe inter- 
ference effect a 23 ' 24 ' 25 ' 26 , e.g. in conjugated monocyclic 
molecules as benzene or annulen o 23 ' 25 . However, for 
the description of transport through a molecule weakly 
coupled to leads, other methods are required. In the 
Coulomb Blockade regime, for example, due to the cru- 
cial role played by the Coulomb interaction in these sys- 
tems it is common to resort to a Pauli rate equation 2 - 7 - 
or to a generalized master equation for the reduced den- 
sity matrix (RDM). For example, in Hettler et al£&, an 
electronic structure calculation was performed in order 
to construct an effective interacting Hamiltonian for the 
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7r orbitals of benzene, and the I-V characteristics of the 
corresponding molecular junction were calculated within 
the rate equation approach. 

In the presence of degenerate states, however, coherences 
of the density matrix influence the dynamics and a mas- 
ter equation approach is appropriat e) 29 ' 30 ' 31 ' 32 ' 33 ' 34 ' 35 ' 36 . 
Such coherences can give rise to precession effects in spin 
transpor t 30 ' 3 or cause interference in a molecular single 
electron transisto r 32 ' 34 ' 35 . In the present work we wish 
to generalize the discussion on interference phenomena in 
a benzene interference- SET presented in^ 5 . to the case in 
which the perfect degeneracy is broken due e.g. to con- 
tact effects or to the applied external bias. To this ex- 
tent the master equation used in^ 5 . will be generalized to 
treat the case of quasi-degenerate states. Conditions for 
the persistence of interference phenomena are identified. 
We observe that the effects of quasi-degenerate states on 
transport have been very recently addressed also in^. 
We treat the transport through the benzene I-SET in 
two different setups, the para and the meta configura- 
tion, depending on the position of the leads with respect 
to the benzene molecule (see Fig. [I]). Similar to2£, we 
start from an interacting Hamiltonian of isolated benzene 
where only the localized p z orbitals are considered and 
the ions are assumed to have the same spatial symmetry 
as the relevant electrons. We calculate the 4 6 = 4096 en- 
ergy eigenstates of the benzene Hamiltonian numerically. 

Subsequently, with the help of group theory, we classify 
the eigenstates according to their different symmetries 
and thus give a group-theoretical explanation to the large 
degeneracies occurring between the electronic states. For 
example, while the six-particles ground state (A\ g sym- 
metry) is non-degenerate, there exist four seven-particle 
ground states due to spin and orbital (i?2« symmetry) 
degeneracy. Fingerprints of these orbital symmetries are 
clearly visible in the strong differences in the stability di- 
agrams obtained by coupling the benzene I-SET to the 
leads in the meta and para configurations. Striking are 
the selective reduction of conductance and the appear- 
ance of regions of interference driven current blocking 
with associated negative differential conductance (NDC) 
when changing from the para to the meta configuration. 

NDC and current blocking for benzene junctions have 
been predicted also in 2 ^, but in the para configuration 
and in presence of an external electromagnetic field. In 
our work NDC occurs despite the absence of an external 
field in the unperturbed setup and with no asymmetry 
in the tunnelling rates. In fact, NDC and current block- 
ing triggered by interference take place any time a SET 
presents an iV-particle non-degenerate state and two de- 
generate AT + 1-particle states such that the ratio between 
the transition amplitudes ji a (i = 1,2, a — L,R) be- 
tween those N- and N + 1-particle states is different for 
tunneling at the left (L) and at the right (R) lead: 

m^TiH. (i) 

Notice that no asymmetry in the tunnelling rates, which 



are proportional to 7i a | 2 , is implied by (JTJ) . This fact 
excludes the interpretation of the physics of the interfer- 
ence SET in terms of standard NDC with asymmetric 
couplings. Due to condition (1) there exist linear com- 
binations of the degenerate N + 1-particle states which 
are coupled to one of the leads but not to the other. The 
state which is decoupled from the right lead represents a 
blocking state for the current flowing L — > R since elec- 
trons can populate this state by tunnelling from the left 
lead but cannot tunnel out towards the right lead. Vicev- 
ersa the state decoupled from the left lead is a blocking 
state for the current R —> L. Typically these two blocking 
states are not orthogonal and thus cannot form together 
a valid basis set. The basis set that diagonalizes the sta- 
tionary density matrix (what we call in the manuscript 
the "physical basis") contains at large positive biases the 
L — > R blocking state and is thus different from the 
physical basis at large negative biases which necessar- 
ily contains the R — ► L blocking state. More generally 
the "physical basis" depends continuously on the bias. 
Thus only a treatment that includes coherences in the 
density matrix can capture the full picture at all biases. 
By neglecting for simplicity the spin degree of freedom, 
the 7-particle ground state of benzene is two times de- 
generate while the 6-particle one is non-degenerate. If 
we choose for the 7-particle states the eigenstates of the 
z-projection of the angular momentum we obtain the re- 
lation: 

Til = 7rR e 4i0 ^) 

72 L J2R 

where 4> is the angle between the left and the right lead. 
Thus in the meta configuration (<j> = 2ir/3) the condi- 
tion is fullfilled while in the para {4> = n) the am- 
plitude ratios are equal. This condition implies that, in 
the para configuration one of the 7-particle states is de- 
coupled from both leads at the same time and can thus 
(in first approximation) be excluded from the dynam- 
ics. In contrast, in the meta configuration, the linear 
combination of uniformly distributed eigenstates of the 
angular momentum creates states with a peculiar inter- 
ference pattern. The position of their nodes allows to 
characterize them as different blocking states. 

This paper is outlined as follows: in Section[TT]we intro- 
duce the model Hamiltonian of the system and present a 
density matrix approach setting up a generalized master 
equation describing the electron dynamics. We give the 
expression for the current in the fully symmetric setup 
(the generalized master equation and current formula for 
the setup under perturbation are given in Appendix [A]) . 
Further we provide a detailed analysis of the symmetry 
characteristics of the molecular eigenstates. 
In Section |Hll we present numerical and analytical results 
of transport calculations for the unperturbed setup. We 
study the occurring interference effects and provide an 
explanation of the phenomena based on symmetry con- 
siderations. 

In Section IIVI we present the results for the perturbed 
setup including a detailed discussion of the transport in 
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this case. We identify in the quasi- degeneracy of the con- 
tributing molecular states the necessary condition for the 
robustness of the interference effects. 
Conclusions and remarks are presented in section fVl 

II. MODEL HAMILTONIAN AND THE 
DENSITY MATRIX APPROACH 

A. Model Hamiltonian 

For the description of the benzene molecule weakly 
coupled to source and drain leads, we adopt the total 
Hamiltonian H = #g en + Hi cads + H T + H^ cn . The 
first term is the interacting Hamiltonian for isolated 
benzene 3 -^^ 

+£/^(n iT -i)(n i; -|) ( 3 ) 

i 

+V^2(n^ + mi - 1) (n i+ it + m+n - 1) , 

i 

where d\ a creates an electron of spin a in the p z orbital 
of carbon i, i = 1, . . . , 6 runs over the six carbon atoms 
of benzene and rii a = d\ a d ia . 

Only the p z orbitals (one per carbon atom) are explic- 
itly taken into account, while the core electrons and the 
nuclei are combined into frozen ions, with the same spa- 
tial symmetry as the relevant electrons. They contribute 
only to the constant terms of the Hamiltonian and en- 
force particle-hole symmetry. Mechanical oscillations are 
neglected and all atoms are considered at their equilib- 
rium position. 

This Hamiltonian for isolated benzene is respecting the 
Dq^ symmetry of the molecule. Since for every site there 
are 4 different possible configurations (|0), |J.), |tl)), 
the Fock space has the dimension 4 6 — 4096, which re- 
quires a numerical treatment. Though the diagonaliza- 
tion of the Hamiltonian is not a numerical challenge, it 
turns out to be of benefit for the physical understand- 
ing of the transport processes to divide -Hbcn into blocks, 
according to the number ./V of p z electrons (from to 
12), the z projection S z of the total spin and the orbital 
symmetries of benzene (see Table HJ. 
The parameters b, U, and V for isolated benzene are 
given in the literature 4 ^ and are chosen to fit optical ex- 
citation spectra. The presence, in the molecular I-SET, 
of metallic electrodes and the dielectric is expected to 
cause a substantial renormalization of U and V— £L. Nev- 
ertheless, we do not expect the main results of this work 
to be affected by this change. We consider the benzene 
molecule weakly coupled to the leads. Thus, to first ap- 
proximation, we assume the symmetry of the isolated 
molecule not to be changed by the screening. Perturba- 
tions due to the lead-molecule contacts reduce the sym- 
metry in the molecular junction. They are included in 



ffben ( see Ec l- (E3) and ([26])) and will be treated in Sec- 
tion [TV] 

The effect of the gate is included as a renormalization of 
the on-site energy £ = £ — eV g (V g is the gate voltage) 
and we conventionally set V g = at the charge neutrality 
point. Source and drain leads are two reservoirs of non- 
interacting electrons: H lcads = Y,aka( e k - ^ a )^ ak(7 c aka , 
where a = L,R stands for the left or right lead and the 
chemical potentials p a of the leads depend on the ap- 
plied bias voltage /il.r = Mo i 2 ■ I n the following we 
will measure the energy starting from the equilibrium 
chemical potential po = 0. The coupling to source and 
drain leads is described by the tunnelling Hamiltonian 

ffT=iE(4/ ato +4A). ( 4 ) 

aka 

where we define d) aa as the creator of an electron in the 
benzene carbon atom which is closest to the lead a. In 
particular d^ := d 4a ,d 5cr respectively in the para and 
meta configuration, while d\^ a := d\ a in both setups. 

B. Dynamics of the reduced density matrix 

Given the high degeneracy of the spectrum, the 
method of choice to treat the dynamics in the weak cou- 
pling is the Liouville equation method already used e.g. 
^32,33 ^ 

In this section we shortly outline how to derive 
the equation of motion for the reduced density matrix 
(RDM) to lowest non-vanishing order in the tunnelling 
Hamiltonian. For more details we refer to 33 ' 34 . 
Starting point is the Liouville equation for the total den- 
sity operator of molecule and leads p in the interac- 
tion picture, treating Ht as a perturbation: ih dp d ^ = 
[iJy, p l {t)] . This equation integrated over time and iter- 
ated to the second order reads 

p l (t)=-~[H l T (t),p l (t )] (5) 

^ f t dt '[H\{t),[HW),p\t')]]- 

Since we are only interested in the transport through the 
molecule, we treat from now on the time evolution of 
the reduced density matrix (RDM) a = Tricadsjp 1 ^)}^, 
which is formally obtained from Eq. (0 by tracing out 
the leads degrees of freedom: a = Trieadsjp 1 }. 
In order to proceed, we make the following standard ap- 
proximations: 

i) The leads are considered as reservoirs of non- 
interacting electrons in thermal equilibrium. Hence we 
can factorize the density matrix as p l (t) = a{t)p s p d — 

c(£)Plcads- 

ii) Since the moecule is weakly coupled to the leads we 
treat the effects of to the lowest non- vanishing order. 

iii) Due to the continuous interaction of the system with 
the leads and at high enough temperature, it is legiti- 
mate to apply the Markov approximation and obtain an 
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equation for & which is local in time (cr(i) instead of cr(t') 
inside the integral) . In particular the Markov approxima- 
tion becomes exact in the stationary limit (i — > oo) we 
will focus on. Since we are interested in the long term 
behavior of the system, we set to — ► — oo in Eq. |(SJ) and 
finally obtain the generalized master equation (GME) 

<j(t) = j^J *"Tr leads {[/4(t), [fl* (t - t"),a(t)pi eada ]}}. 

(6) 

The reduced density operator a is defined on the Fock 
space of benzene, yet we can neglect coherences between 
states with different particle number since decoupled 
from the dynamics of the populations. For simplicity, 
we continue here the derivation of the GME only for 
the symmetric case with exact orbital degeneracy, i.e., 
neglecting (the perturbed case is presented in Ap- 

pendix [2J . 

iv) Further we also neglect coherences between states 
with different energy (secular approximation). They are 
irrelevant due to their fast fluctuation compared to the 
dynamics of the system triggered by the tunnelling cou- 



pling. 

Under these considerations, it is convenient to ex- 
press the GME in terms of the reduced density operator 
t NE = VnetVne, where T NE := Eir \NE£t)(NE£t\ 
is the projection operator on the subspace of N parti- 
cles and energy E. The sum runs over the orbital and 
spin quantum numbers £ and r, respectively. The orbital 
quantum number £ distinguishes between orbitally degen- 
erate states. The exact meaning of £ will be illustrated in 
the next section. In Appendix lAl we derive a GME that 
retains coherences also between quasi-degenerate states. 
That approach treats with special care the small asym- 
metries introduced in the molecule by the coupling to 
the leads. In fact it interpolates between the degenerate 
case treated here and the fully non-degenerate case in 
which the GME reduces to a master equation for pop- 
ulations only. Equation ^ can be further manipulated 
by projection into the subspace of TV-particle and energy 
E. Since we assume the density matrix to be factorized 
and the leads to be in thermal equilibrium, also the traces 
over the leads degree of freedom can be easily performed. 
Eventually, the GME for the degenerate case reads 



a NE = 



fa ( H bcn - E ) - — Pa (#bcn ~ E ) 



fa ( E ~ H bcn) Pa{E ~ H^ D 



d aT a NE + H.c. 



(7) 



E r ^NE \dtM(E E>) a N - 1E 'd aT + d ar f-(E> - E) CT N+1E '4 T W, 

arE' I J 



where Tl.r = I^l.rI^l.r are the bare transfer rates 
with the constant densities of states of the leads 2?l,r- 
Terms describing sequential tunnelling from and to the 
lead a are proportional to the Fermi functions f^(x) := 
f(x — fi a ) an d fa( x ) : = 1 ~ /a"( x )i respectively. Still 
in the sequential tunnelling limit, but only in the equa- 
tions for the coherences, one finds also the energy non- 
conserving terms, proportional to the function p a (x) — 



-Retp 



i/3 , 



Ma) 



, where if) is the digamma func- 
tion. Both the Fermi functions and the digamma func- 
tion result from the trace over the leads degrees of 
freedom22^£. 

A closer analysis of the master equation allows also to for- 
mulate an expression for the current operator. We start 
from the definition of the time derivative of the charge 
on benzene: 



dt 



(Q) = Tr 



{Na} 



= (Il+Ik) 



(8) 



where Q = ^2 iT {d\ T d iT — 6) is the operator of the charge 
on benzene, N is the particle number operator and 7l,r 
are the current operators at the left(right) contact. Con- 
ventionally, in the definition of 7l,r we assume the cur- 
rent to be positive when it is increasing the charge on the 
molecule. Thus, in the stationary limit, ( Zl +/r ) is zero. 
We write this expression in the basis of the subspaces of 
N particles and energy E: 



(Il + Ir)=J2 Tr {^nk<tPne} = E Tr {^ NE } ■ 

NE NE 

(9) 

Further we insert {7} in ([9]) and take advantage of the 
cyclic properties of the trace to find : 



5 



i J l + In > = E E - [fi( H ln E)d ar a NE d aT + f-(E Hl n )d aT a^dt T ] 

NE ar { 



(10) 



5>ne [/+(£ - S')4^ N - 1E '^ + fa (E' E)d aT a^'di T 

E' 

I 



Notice that the energy non-conserving contributions drop 
from the expression of the current. Still they contribute 
to the average current via the density matrix. Since E 
and E' are dummy variables, we can switch them in the 
summands containing E' . Applying the relation: 



NE' 



Tr{VKwg(E')} = Tr{g(H° en )}, 



where g{E') is a generic function, we substitute E' with 
-^ben m Eq. pop . Further we can conveniently rearrange 
the sum over N, arriving at the expression for the cur- 
rent: 



;/ L +/ R > = EE r « Tr 

NE ar 

dt T a* E d\-Nf+(H° cn -E) + (N + l)f+(Hl n -E) 



-da NE dt 



-Nf-(E-Hg en ) + (N-l)f-(E- HU\ | • 

(11) 



This relation can be further simplified in order to iden- 
tify the current operators. The one corresponding to the 
left contact is e.g. 



NEt 



/l = T L E ^ne [d hT ti(H° cn - E)d{ T + 



(12) 



Vne- 



With this relation we can calculate the stationary current 
as the average (II) = Tr{cr sta tiL} = — (Ir), with cr sta t 
as the stationary density operator. The expression of 
the current operator for the perturbed system is given in 
Appendix [A] 



C. Symmetry of the benzene eigenstates 

In this section, we will review the symmetry charac- 
teristics of the eigenstates of the interacting Hamiltonian 
of benzene, focusing on the symmetry operations a v and 
C n which have a major impact on the electronic trans- 
port through the molecular I-SET. Benzene belongs to 
the D^h point group. Depending on their behavior un- 
der symmetry operations, one can classify the molecular 



orbitals by their belonging to a certain irreducible repre- 
sentation of the point group. 



N #T #1 # states 


# states with a 
certain symmetry 


6 6 1 


1 B lu 


5 1 36 


2A 2s 
2x6 E2 9 
4 B lu 
2 B 2 u 
2x6 E lu 


4 2 225 


16 A lg 
20 A 2g 

2x36 E 2g 
22 Biu 

17 B 2u 
2x39 E lu 


3 3 400 


38 A lg 
30 A 2g 

2x66 E 2g 
38 Biu 
30 B 2u 

2x66 Eiu 


2 4 225 

1 5 36 
6 1 





TABLE I: Overview of the 6 particle states of benzene, sorted 
by S z and symmetry. Orbitals with A- and B-type of sym- 
metry show no degeneracy, while _B-type orbitals are doubly 
degenerate. 

Table U shows an overview of the states of the neutral 
molecule (the 6 particle states) sorted by S z and symme- 
tries. The eigenstates of the interacting benzene molecule 
have either A-, B- or E-type symmetries. While orbitals 
having A or B symmetries can only be spin degenerate, 
states with an E symmetry show an additional twofold 
orbital degeneracy, essential for the explanation of the 
transport features occurring in the meta configuration. 

Transport at low bias is described in terms of transi- 
tions between ground states with diferent particle num- 
ber. Table |n] shows the symmetries of the ground states 
(and of some first excited states) of interacting ben- 
zene for all possible particle numbers. Ground state 
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transitions occur both between orbitally non-degenerate 
states (with A and B symmetry), as well as between or- 
bitally degenerate and non-degenerate states (E- to A- 
type states). 

The interacting benzene Hamiltonian commutes with 
all the symmetry operations of the D§h point group, thus 
it has a set of common eigenvectors with each operation. 
The element of Dq^ of special interest for the para config- 
uration is tr v , i.e., the reflection about the plane through 
the contact atoms and perpendicular to the molecular 
plane. The molecular orbitals with A and B symmetry 
are eigenstates of a v with eigenvalue ±1, i.e., they are 
either symmetric or antisymmetric with respect to the 
<7 V operation. The behavior of the E— type orbitals un- 
der a v is basis dependent, yet one can always choose a 
basis in which one orbital is symmetric and the other one 
antisymmetric. 

Let us now consider the generic transition amplitude 
(N\d aT \N + 1), where d aT destroys an electron of spin r 
on the contact atom closest to the a lead. It is useful to 
rewrite this amplitude in the form 

(N\d aT \N+l) = (N\<T$a v d aT al* v \N + l), (13) 

where we have used the property crjer v = 1. Since in the 
para configuration both contact atoms lie in the mirror 
plane <r v , it follows a w d a a\ — d a . If the participating 
states are both symmetric under cr v , Eq. (|13[) becomes 

(N, sym\ald aT a v \N + 1, sym) = 

= (N,sym\d aT \N + l,sym). (14) 

and analogously in the case that both states are antisym- 
metric. For states with different symmetry it is 

(N, sym|ci aT |JV + 1, antisym) = 

= — (N, sym\d aT \N + 1, antisym) =0. (15) 

In other terms, there is a selection rule that forbids 
transitions between symmetric and antisymmetric states. 
Further, since the ground state of the neutral molecule is 
symmetric, for the transport calculations in the para con- 
figuration we select the effective Hilbert space containing 
only states symmetric with respect to a v . Correspond- 
ingly, when referring to the N particle ground state we 
mean the energetically lowest symmetric state. For ex- 
ample in the case of 4 and 8 particle states it is the first 
excited state to be the effective ground state. In the para 
configuration also the orbital degeneracy of the E— type 
states is effectively cancelled due to the selection of the 
symmetric orbital (see Table ITT|) . 

Small violations of this selection rule, due e.g. to 
molecular vibrations or coupling to an electromag- 
netic bath, result in the weak connection of different 
metastable electronic subspaces. We suggest this mech- 
anism as a possible explanation for the switching and 
hysteretic behaviour reported in various molecular junc- 
tions. This effect is not addressed in this work. 

For a simpler analysis of the different transport char- 
acteristics it is useful to introduce a unified geometrical 



description of the two configurations. In both cases, one 
lead is rotated by an angle <f> with respect to the posi- 
tion of the other lead. Hence we can write the creator 
of an electron in the right contact atom d\^ T in terms of 
the creation operator of the left contact atom and the 
rotation operator: 

4r = K d l^ (16) 

where TZ,p is the rotation operator for the anticlockwise 
rotation of an angle <p around the axis perpendicular to 
the molecular plane and piercing the center of the ben- 
zene ring; <p = 7T for the para and <f> = (2ir/3) for the 
meta configuration. 

N Degeneracy Energy[eV] Symmetry Symmetry behavior 
(at f = 0) under ov 

1 A lg sym 

1 2 -22 A 2u sym 

2 1 -42.25 Aig sym 

3 4 -57.42 Ei g 2 sym, [2 antisym] 

4 [3] [-68.87] [A 2g ] [antisym] 
2 -68.37 E 2g 1 sym, [1 antisym] 

5 4 -76.675 E lg 2 sym, [2 antisym] 

6 1 -81.725 Aig sym 

7 4 -76.675 E 2u 2 sym, [2 antisym] 

8 [3] [-68.87] [Aug] [antisym] 
2 -68.37 E 2g 1 sym, [1 antisym] 

9 4 -57.42 E 2u 2 sym, [2 antisym] 

10 1 -42.25 A lg sym 

11 2 -22 B 2g sym 

12 1 A lg sym 

TABLE II: Degeneracy, energy and symmetry of the ground 
states of the isolated benzene molecule for different particle 
numbers. We choose the on-site and inter-site Coulomb in- 
teractions to be U — 10 eV, V = 6eV, and the hopping to 
be b — — 2.5eV. Notice, however, that screening effects from 
the leads and the dielectric are expected to renormalize the 
energy of the benzene many-body states. 

The energy eigenstates of the interacting Hamiltonian 
of benzene can be classified also in terms their quasi- 
angular momentum. In particular, the eigenstates of the 
z-projection of the quasi angular momentum are the ones 
that diagonalizc all operators TZ^ the with angles multi- 
ples of 7r/3. The corresponding eigenvalues are phase fac- 
tors e~ %1 ^ where M, the quasi-angular momentum of the 
state, is an integer multiple of h. The discrete rotation 
operator of an angle (f> = n (C2 symmetry operation), is 
the one relevant for the para configuration. All orbitals 
are eigenstates of the C2 rotation with the eigenvalue ±1. 

The relevant rotation operator for the meta configu- 
ration correspond to an angle (f> = 2ir/3 (C3 symmetry 
operation). Orbitals with an A or B symmetry are eigen- 
states of this operator with the eigenvalue +1 (angular 
momentum £ = or £ — 3). Hence we can already pre- 
dict that there will be no difference based on rotational 
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Alg 


sym 


2 


-22 


A 2u 


sym 


1 


-42.25 


Alg 
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symmetry between the para and the meta configuration 
for transitions between states involving A- and B-type 
symmetries. Orbitals with E symmetry however behave 
quite differently under the C3 operation. They are the 
pairs of states of angular momenta I — ±1 or I = ±2. 
The diagonal form of the rotation operator on the two- 
fold degenerate subspace of .©-symmetry reads: 



C, 








(17) 



For the two-fold orbitally degenerate 7-particle ground 
states \£\ = 2. This analysis in terms of the quasi- angular 
momentum makes easier the calculation of the fundamen- 
tal interference condition @ given in the introduction. 
In fact the following relation holds between the transition 
amplitudes of the 6 and 7 particle ground states: 



lm = (V r l<4rl 6 s) 



= {T g ir\Tl\ci r n^ |6 ff > = e~ 



(18) 



and ([2|) follows directly. 
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FIG. 2: (color online) Stability diagram for the benzene I- 
SET contacted in the para (above) and meta (below) con- 
figuration. Dot-dashed lines highlight the conductance cuts 
presented in Fig. [3] the dashed lines the regions corresponding 
to the current traces presented in Fig.|4]and Fig. [6] the dotted 
line the region corresponding to the current trace presented in 
Fig. [5] Fhe parameters used are U = 4\b\, V = 2.4|fe|, k B T = 
0.04|6|, fir L = hV R = 10 3 |fo . 



III. TRANSPORT CALCULATIONS: FULLY 
SYMMETRIC SETUP 

With the knowledge of the eigenstates and eigenvalues 
of the Hamiltonian for the isolated molecule, we imple- 
ment Eq. ([7]) and look for a stationary solution. The 
symmetries of the eigenstates are reflected in the transi- 
tion amplitudes contained in the GME. We find numeri- 
cally its stationary solution and calculate the current and 
the differential conductance of the device. In Fig. [5] we 
present the stability diagram for the benzene I-SET con- 
tacted in the para (upper panel) and meta position (lower 
panel). Bright ground state transition lines delimit di- 
amonds of zero differential conductance typical for the 
Coulomb blockade regime, while a rich pattern of satellite 
lines represents the transitions between excited states. 
Though several differences can be noticed, most striking 
are the suppression of the linear conductance, the appear- 
ance of negative differential conductance (NDC) and the 
strong suppression of the current at the right (left) border 
of the 7 (5) particle diamond when passing from the para 
to the meta configuration. All these features are differ- 
ent manifestations of the interference between orbitally 
degenerate states and ultimately reveal the specific sym- 
metry of benzene. 



A. Linear conductance 

We study the linear transport regime both numerically 
and analytically. For the analytical calculation of the 
conductance we consider the low temperature limit where 



only ground states with N and N + 1 particles have con- 
siderable occupation probabilities, with N fixed by the 
gate voltage. Therefore only transitions between these 
states are relevant and we can treat just the terms of Q 
with N and N + 1 particles and the ground state ener- 
gies _E gi N and 2? g) N+i, respectively. A closer look to (O 
reveals that the spin coherences are decoupled from the 
other elements of the density matrix. Thus we can set 
them to zero, and write ([7]) in a block diagonal form in 
the basis of the ground states of N and N+l particles. 
Additionally, since the total Hamiltonian H is symmet- 
ric in spin, the blocks of the GME with the same particle 
but different spin quantum number r must be identical. 
Finally, since around the resonance the only populated 
states are the N and N+l particle states, the conserva- 
tion of probability implies that: 

! = £<£. + £<££, (19) 

n tci 

where cr„„ is the population of the N-particle ground 
state and n contains the orbital and spin quantum num- 
bers. With all these observations we can reduce |7J to a 
much smaller set of coupled differential equations, that 
can be treated analytically. The stationary solution of 
this set of equations can be derived more easily by ne- 
glecting the energy non-conserving terms in (J7J) . These 
are contained in the elements of the GME describing the 
dynamics of the coherences between orbitally degenerate 
states. With this simplification we derive an analyti- 
cal formula for the conductance close to the resonance 
between N and iV + 1 particle states as the first order 
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coefficient of the Taylor series of the current in the bias: 

r L r R 

-An,n+i x 

(20) 



G 



N.N+l 



(AE) =2e 2 



r R "'"^ 
S N S N+1 f'(AE) 



(5 N+ i - S N )f(AE) + 5^ 



where A£J = S g ,N — -Eg,N+i + eVJ, is the energy difference 
between the benzene ground states with N and N + 1 
electrons diminished by a term linear in the gate voltage. 
Interference effects are contained in the overlap factor 

Aw N+l : 



A 



E (JV,n|d LT |iV+l,mXJV+l,m|4,Tl^.»> 



N,N+1 



SnSn+i E (N,n\d aT \N+l,m) 



where n and m label the SN-fold and iSw+i-fold degen- 
erate ground states with N and N + 1 particles, respec- 
tively. In order to make the interference effects more 
visible we remind that c? Rt = IZld^IZi, with <j) = n for 
the para while 4> — 27r/3 for the meta configuration. Due 
to the behaviour of all eigenstates of H^ cn under discrete 
rotation operators with angles multiples of 7r/3, we can 
rewrite the overlap factor: 



transitions. The ground states with N — 10, 11, 12 par- 
ticles have A— or B— type symmetries, they are there- 
fore orbitally non-degenerate, no interference can occur 
and thus the transitions are invariant under configuration 
change. For every other transition we see a noticeable 
difference between the results of the two configurations 
(Fig. [3]). In all these transitions one of the participat- 
ing states is orbitally degenerate. First we notice that 
the linear conductance peaks for the 7 <-> 8 and 8^9 
transitions in the para configuration are shifted with re- 
spect to the corresponding peaks in the meta configu- 
ration. The selection of an effective symmetric Hilbert 
space associated to the para configuration reduces the 
total degeneracy by cancelling the orbital degeneracy. In 
addition, the ground state energy of the 4 and 8 particle 
states is different in the two configurations, since in the 
para configuration the effective ground state is in reality 
the first excited state. The degeneracies 5n,-Sn+i of the 
participating states as well as the ground state energy 
are both entering the degeneracy term of Eq. (|20|) 



A 



/'(A£) 



(S^-S^f^ + SK' 
and determine the shift of the conductance peaks. 



(22) 



A 



E \(N,n\d LT \N + l,m)\< 



N,N+1 



S N S N+1 2 E (N,n\d LT \N+l,m) 



r, (21) 



where (j) nrn encloses the phase factors coming from the 
rotation of the states \N, n) and \N + 1, m). 

The energy non-conserving terms neglected in (|20p in- 
fluence only the dynamics of the coherences between or- 
bitally degenerate states. Thus, Eq. (|2Tj|) provides an 
exact description of transport for the para configuration, 
where orbital degeneracy is cancelled. Even if Eq. (j20|) 
captures the essential mechanism responsible for the con- 
ductance suppression, we have derived an exact ana- 
lytical formula also for the meta configuration and we 
present it in Appendix [Bl 

In Fig. [3] we present an overview of the results of both 
the para and the meta configuration. A direct compari- 
son of the conductance (including energy non-conserving 
terms) in the two configurations is displayed in the up- 
per panel. The lower panel illustrates the effect of the 
energy non-conserving terms on the conductance in the 
meta configuration. The number of p z electrons on the 
molecule and the symmetry of the lowest energy states 
corresponding to the conductance valleys are reported. 
The symmetries displayed in the upper panel belong to 
the (effective) ground states in the para configuration, 
the corresponding symmetries for the meta configuration 
are shown in the lower panel. 

Fig. [3J shows that the results for the para and the meta 
configuration coincide for the 10 <-> 11 and 11 <-> 12 
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FIG. 3: (color online) Conductance of the benzene I-SET as a 
function of the gate voltage. Clearly visible are the peaks cor- 
responding to the transitions between ground states with TV 
and N + 1 particles. In the low conductance valleys the state 
of the system has a definite number of particles and symme- 
try as indicated in the upper panel for the para, in the lower 
for the meta configuration. Selective conductance suppres- 
sion when changing from the meta to the para configuration 
is observed. 

Yet, the most striking effect regarding transitions with 
orbitally degenerate states participating is the system- 
atic suppression of the linear conductance when changing 
from the para to the meta configuration. The suppres- 
sion is appreciable despite the conductance enhancement 
due to the energy non-conserving terms (see Fig. [31 lower 
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panel). Thus, we will for simplicity discard them in the 
following discussion. 

The conductance suppression is determined by the 
combination of two effects: the reduction to the symmet- 
ric Hilbert space in the para configuration and the inter- 
ference effects between degenerate orbitals in the meta 
configuration. The reduction to the symmetric Hilbert 
space implies also a lower number of conducting chan- 
nels (see Table IIII[) . One would expect a suppression of 
transport in the para configuration. As we can see from 
Table IIIII on the example of the 6^7 transition peak, 
A max is higher in the para configuration but not enough 
to fully explain the difference between the two configu- 
rations. 



B. Negative differential conductance (NDC) and 
current blocking 

Interference effects between orbitally degenerate states 
are also affecting non-linear transport and producing in 
the meta configuration current blocking and thus NDC 
at the border of the 6 particle state diamond (Fig. [2]). 
The upper panel of Fig. [4] shows the current through 
the benzene I-SET contacted in the meta configuration 
as a function of the bias voltage. The current is given 
for parameters corresponding to the white dashed line of 
Fig. [21 In this region only the 6 and 7 particle ground 
states are populated. 





# Channels 


Overlap factor 


Degeneracy term 




SnSn+i 


A 


A max [1/fcsT] 


PARA 


2 


C 


0,17 


META 


4 


\c 


0,11 



TABLE III: Number of channels participating to transport, 
overlap factor and resonance value of the degeneracy term in 
the para and the meta configuration for the 6^7 transition 
peak. It is G = | (6 g Ic^lt- 1 7 g Zx) | 2 , where r and £ are the spin 
and the quasi angular momentum quantum numbers, respec- 
tively. The values of A max are given for ksT — 0.04|&| 



The second effect determining transport is the interfer- 
ence between the -E-type states, which is accounted for 
in the overlap factor A. The overlap factor is basis in- 
dependent, thus we can write the transition probabilities 
for the 6^7 transition as |(6g|d Lr |7 g £T}| 2 = C, where 
r and I are the spin and the quasi-angular momentum 
quantum number, respectively. The transition probabil- 
ities have the same value, since all four 7 particle states 
are in this basis equivalent (see Appendix [C|- Under the 
C*2 rotation the symmetric 7 particle ground state does 
not acquire any phase factor. Under the C3 rotation how- 
ever, the two orbitally degenerate states acquire different 
phase factors, namely e~ l and e~~ l , respectively. Thus 
the overlap factors A for the 6 «-> 7 transition are: 



A 



pare 



1 |4C| 2 = C, 



A n 



2-8C 
1 

4-8C 



2Ce + ~ l + 2Ce~- 



= lc. 



The linear conductance is determined by the product be- 
tween the number of conducting channels, the overlap 
factor and the degeneracy term. Yet, it is the destruc- 
tive interference between degenerate E-type orbitals, ac- 
counted for in the overlap factor A, that gives the major 
contribution to the strong suppression of the conductance 
in the meta configuration. 




Resonance 
7g Jg_ 

^6, / 



Blocking 
7g 7g 



FIG. 4: (color online) Upper panel - Current through the ben- 
zene I-SET in the meta configuration calculated at bias and 
gate voltage conditions indicated by the dashed line of Fig. [2] 
A pronounced NDC with current blocking is visible. Middle 
panels - Transition probabilities between the 6 particle and 
each of the two 7 particle ground states for bias voltage val- 
ues labelled a — c in the upper panel. The transition to a 
blocking state is visible in the upper (lower) part of the c (a) 
panels. Lower panels - Sketch of the energetics for the 6^7 
transition in the meta configuration at bias voltages corre- 
sponding to the resonance current peak and current blocking 
as indicated in the upper panel of this figure. 
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At low bias the 6 particle state is mainly occupied. As 
the bias is raised, transitions 6 «-> 7 occur and current 
flows. Above a certain bias threshold a blocking state is 
populated and the current drops. For the understand- 
ing of this non- linear current characteristics, we have to 
take into account energy conservation, the Pauli exclu- 
sion principle and the interference between participating 
states. For the visualization of the interference effects, we 
introduce the transition probability (averaged over the z 
coordinate and the spin a): 

P(x,y;n, r) = lim £ -L / L/2 dz|(7 g nr|^(r)|6 g )| 2 

for the physical 7 particle basis, i.e., the 7 particle ba- 
sis that diagonalizes the stationary density matrix at a 
fixed bias. Here t is the spin quantum number, n = 1,2 
labels the two states of the physical basis which are lin- 
ear combinations of the orbitally degenerate states \7 g £r) 
and can be interpreted as conduction channels. Each of 
the central panels of Fig. [4] are surface plots of (|23|) at 
the different bias voltages a-c. The 7 particle ground 
states can interfere and thus generate nodes in the tran- 
sition probability at the contact atom close to one or the 
other lead, but, in the meta configuration, never at both 
contact atoms at the same time. 

Energetic considerations are illustrated in the lower 
panels of Fig. U for two key points of the current curve 
at positive biases. The left panel corresponds to the res- 
onance peak of the current. Due to energy conservation, 
electrons can enter the molecule only from the left lead. 
On the contrary the exit is allowed at both leads. The 
current is suppressed when transitions occur to a state 
which cannot be depopulated (a blocking state). Since, 
energetically, transmissions to the 6 particle state are al- 
lowed at both leads, each 7 particle state can always be 
depopulated and no blocking occurs. 

The current blocking scenario is depicted in the lower 
right panel of Fig. |H For large positive bias the transition 
from a 7 particle ground state to the 6 particle ground 
state is energetically forbidden at the left lead. Thus, 
for example, the c panel in Fig. [4] visualizes the current 
blocking situation yielding NDC: while for both channels 
there is a non-vanishing transition probability from the 
source lead to the molecule, for the upper channel a node 
prevents an electron from exiting to the drain lead. In 
the long time limit the blocking state gets fully populated 
while the non-blocking state is empty. At large negative 
bias the blocking scenario is depicted in the panel a that 
shows the left-right symmetry obtained by a reflection 
through a plane perpendicular to the molecule and pass- 
ing through the carbon atoms atoms 6 and 3. 
We remark that only a description that retains coher- 
ences between the degenerate 7 particle ground states 
correctly captures NDC at both positive and negative 
bias. 

In contrast to the 6^7 transition, one does not ob- 
serve NDC at the border of the 7 particle Coulomb di- 



amond, but rather a strong suppression of the current. 
The upper panel of Fig. [5] shows the current through the 
benzene I-SET contacted in the meta configuration as a 
function of the bias voltage corresponding to the white 
dotted line of Fig. [2] The middle panels show the transi- 
tion probabilities between each of the 7 particle and the 
6 particle ground state. 




Resonance 




FIG. 5: (color online) Upper panel - Current through the 
benzene I-SET in the meta configuration calculated at bias 
and gate voltage conditions indicated by the dotted line of Fig. 
[2] No NDC is visible. Middle panels - Transition probabilities 
between each of the 7 particle and the 6 particle ground state 
for bias voltage values labelled a — cm the upper panel. Lower 
panel - Sketch of the energetics for the 7^6 transition in 
the meta configuration at bias voltage corresponding to the 
expected resonance peak, (compare to Fig. [4}. 

The lower panel of Fig. [5] shows a sketch of the ener- 
getics at positive bias corresponding to the "expected" 
resonance peak. Here electrons can enter the molecular 
dot at both leads, while the exit is energetically forbidden 
at the left lead. Thus, if the system is in the 7 particle 
state which is blocking the right lead, this state cannot 
be depopulated, becoming the blocking state. 
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On the other hand, transitions from the 6 particle ground 
state to both 7 particle ground states are equally proba- 
ble. Thus the blocking state will surely be populated at 
some time. The upper plot of the b panel in Fig. \5\ shows 
the transition probability to the blocking state that ac- 
cepts electrons from the source lead but cannot release 
electrons to the drain. 

As just proved, in this case the current blocking situa- 
tion occurs already at the resonance bias voltage. For a 
higher positive bias, the transition probability from the 
blocking state at the drain lead increases and current can 
flow. This effect, though, can be captured only by taking 
into account also the energy non-conserving terms in 10 . 

In the para configuration, the current as a function of 
the bias voltage is shown in Fig. [6] The current is given 
for parameters corresponding to the white dashed line of 
Fig. [2] In this case, no interference effects are visible. We 
see instead the typical step-like behavior of the current 
in the Coulomb blockade regime. 
The panels on the right are the surface plots of 

1 f L/2 

L^oo^ ZLJ _ L / 2 

The upper plot shows the transition probability to the 
symmetric 7 particle state, the lower to the antisym- 
metric. Remember that in the para configuration only 
the symmetric states contribute to transport. Evidently 
the symmetric state is in the para configuration non- 
blocking. Additionally, since the coherences between or- 
bitally degenerate states and therefore the energy non- 
conserving terms do not play any role in the transport, 
the physical basis states are not bias dependent. Thus 
in the para configuration there are always non-blocking 
states populated and no NDC can occur. 




FIG. 6: (color online) Left panel - Current through the ben- 
zene I-SET in the para configuration calculated at bias and 
gate voltage conditions indicated by the dashed line of Fig. [2] 
No interference effects are visible. Right panels - Transition 
probabilities between the 6 particle and the symmetric and 
antisymmetric 7 particle ground states. 



IV. REDUCED SYMMETRY 

In this section we study the effect of reduced symme- 
try on the results presented previously. We generalize 
the model Hamiltonian by taking into account the per- 
turbations on the molecule due to the contacts and the 
bias voltage. The contact between molecule and leads is 
provided by different anchor groups. These linkers are 
coupled to the contact carbon atoms over a a bond thus 
replacing the corresponding benzene hydrogen atoms. 
Due to the orthogonality of 7r and a orbitals, the anchor 
groups affect in first approximation only the a orbitals 
of benzene. In particular the different electron affinity of 
the atoms in the linkers imply a redistribution of the den- 
sity of a electrons. Assuming that transport is carried by 
7r electrons only, we model the effect of this redistribu- 
tion as a change in the on-site energy for the p z orbitals 
of the contact carbon atoms: 

H Ln : = ^contact = ^ d aa d a *i a= L,R (25) 

Ct<7 

where R — 4, 5, respectively, in the para and meta con- 
figuration, L = 1 in both setups. 

We also study the effect of an external bias on the ben- 
zene I-SET. In particular we release the strict condition 
of potential drop all concentrated at the lead-molecule 
interface. Nevertheless, due to the weak coupling of 
the molecule to the leads, we assume that only a frac- 
tion of the bias potential drops across the molecule. For 
this residual potential we take the linear approximation 
Pb (r) = — -^(r • fsd/ao); where we choose the center of 
the molecule as the origin and r S( j is the unity vector di- 
rected along the source to drain direction, ao = 1.43 A 
is the bond length between two carbon atoms in ben- 
zene, a is the coefficient determining the intensity of the 
potential drop over the molecule. Since the p z orbitals 
are strongly localized, we can assume that this potential 
will not affect the inter-site hopping, but only the on-site 
term of the Hamiltonian: 

#L„:=#bias = e $>AU CT (26) 

i<j 

with £ bi = Jdrp z (r- Ri)V h {r)p z (r - Ro- 
under the influence of the contacts or the bias poten- 
tial, the symmetry of the molecule changes. Table IIVI 
shows the point groups to which the molecule belongs in 
the perturbed setup. This point groups have only A- and 
B-type reducible representations. Thus the correspond- 
ing molecular orbitals do not exhibit orbital degeneracy. 

No interference effects influence the transport in the 
para configuration. Thus we do not expect its transport 
characteristics to be qualitatively modified by the new 
set up with the corresponding loss of degeneracies. 

In the meta configuration on the other hand, interfer- 
ences between orbitally degenerate states play a crucial 
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role in the explanation of the occurring transport fea- 
tures. Naively one would therefore expect that neither 
conductance suppression nor NDC and current blocking 
occur in a benzene I-SET with reduced symmetry. Yet we 





PARA 


META 


Contact perturb. 




C2v 


Bias perturb. 


C2v 





TABbE IV: Point groups to which the molecule belongs under 
the influence of the contacts and the external bias potential. 

find that, under certain conditions, the mentioned trans- 
port features are robust under the lowered symmetry. 

The perturbations due to the contacts and the bias 
lead to an expected level splitting of the former orbitally 
degenerate states. Very different current-voltage charac- 
teristics are obtained depending of the relation between 
the energy splitting SE and other two important energy 
scales of the system: the tunnelling rate T and the tem- 
perature T. In particular, when SE <C I « T, interfer- 
ence phenomena persist. In contrast when Y < SE <C T 
interference phenomena disappear, despite the fact that, 
due to temperature broadening, the two states still can 
not be resolved. In this regime, due to the asymmetry 
in the tunnelling rates introduced by the perturbation, 
standard NDC phenomena, see Fig. [5J occur. 

In the absence of perfect degeneracy, we abandon 
the strict secular approximation scheme that would dis- 
card the coherences in the density matrix between states 
with different energies. We adopt instead a softer ap- 
proximation by retaining also coherences between quasi- 
degenerate states. Since they have Bohr frequencies com- 
parable to the tunnelling rate, they influence the sta- 
tionary density matrix. Formulas for the GME and the 
current taking into account these coherence terms are 
presented in Appendix [Al 

Fig. [7] shows from left to right closeup views of the 
stability diagram for the setup under the influence of 
increasing contact perturbation around the 6 <-> 7 res- 
onance. The orbital degeneracy of the 7 particle states 
is lifted and the transport behavior for the 6^7 tran- 
sition depends on the energy difference between the for- 
merly degenerate 7 particle ground states. In panel a the 
energy difference is so small that the states are quasi- 
degenerate: SE <C HT <C UbT. As expected, we recover 
NDC at the border of the 6 particle diamond and current 
suppression at the border of the 7 particle diamond, like 
in the unperturbed setup. 

Higher on-site energy-shifts correspond to a larger level 
spacing. Panel b displays the situation in which the latter 
is of the order of the level broadening, but still smaller 
than the thermal energy (SE ~ HT <C k^T): no interfer- 
ence causing NDC and current blocking can occur. Yet, 
due to thermal broadening, we cannot resolve the two 7 
particle states. 

Eventually, panel c presents the stability diagram for the 
case SE > k^T > HT: the level spacing between the 7 




FIG. 7: (color online) Closeup views of the stability diagram 
around the 6^7 resonance for the system under contact 
perturbation. The perturbation strength grows from left to 
right The parameter that describes the contact effect assumes 
the values £ c = 0.15r, 2F, 15T from left to right respectively 
and T = 10T . 

particle ground and first excited state is now bigger than 
the thermal energy, thus the two transition lines corre- 
sponding to these states are clearly visible at the border 
of the 6 particle stability diamond. 




FIG. 8: (color online) Closeup views of the stability diagram 
around the 6 <-> 7 resonance for the system under the effect 
of the bias potential, displayed for different strengths of the 
electrostatic potential drop over the molecule. The parameter 
that describe the strength of the electrostatic drop overthe 
molecule assumes the values a = 25, 12, 0.6 from left to right 
respectively. 

Fig. [8J shows closeup views of the stability diagram for 
the setup under the influence of the bias perturbation at 
the border of the 6 and 7 particle diamonds. The same 
region is plotted for different strengths of the external 
potential over the molecule. In contrast to the contact 
perturbation, the amount of level splitting of the former 
degenerate states is here bias dependent. This fact im- 
poses a bias window of interference visibility. The bias 
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must be small enough, for the 7 particle states to be 
quasi-degenerate and at the same time bigger than the 
thermal energy, so that the occurring NDC is not ob- 
scured by the thermally broadened conductance peak. A 
strong electrostatic potential perturbation closes the bias 
window and no interference effect can be detected. 
Panel a of Fig.[5]represents the weak perturbation regime 
with no qualitative differences with the unperturbed case. 
The typical fingerprints of interference (NDC at the bor- 
der of the 6 particle diamond and current blocking for 
the 7^6 transition) are still visible for intermediate 
perturbation strength (panel b) but this time only in a 
limited bias window. Due to the perturbation strength, 
at some point in the bias, the level splitting is so big that 
the quasi-degeneracy is lifted and the interference effects 
destroyed. In panel c the quasi-degeneracy is lifted in 
the entire bias range. There is NDC at the border of the 
6 particle diamond, but is not accompanied by current 
blocking as proved by the excitation line at the border of 
the 7 particle diamond (see arrow): no interference oc- 
curs. The NDC is here associated to the sudden opening 
of a slow current channel, the one involving the 6 particle 
ground state and the 7 particle (non-degenerate) excited 
state (standard NDC). 




p 



FIG. 9: (color online) Combination of the bias and contact 
perturbations. Left panel - Energy levels of the 7 particle 
ground and first excited state as functions of the bias voltage. 
Right panel - Stability diagram around the 6^7 resonance. 
The perturbation parameters are in this case £ c = 2T and 
a = 12. 

Fig. prefers to the setup under both the bias and con- 
tact perturbations. The left panel shows the energy of the 
lowest 7 particle states as a function of the bias. In the 
right panel we present the stability diagram around the 
6 <-> 7 resonance. NDC and current blocking are clearly 
visible only in the bias region where, due to the combina- 
tion of bias and contact perturbation, the two seven par- 
ticle states return quasi-degenerate. Also the fine struc- 
ture in the NDC region is understandable in terms of in- 
terference if in the condition of quasi-degeneracy we take 
into account the renormalization of the level splitting due 



to the energy non-conserving terms. 

Interference effects predicted for the unperturbed ben- 
zene I-SET are robust against various sources of sym- 
metry breaking. Quasi-degeneracy, 6E <C HT -C fceT, 
is the necessary condition required for the detection of 
the interference in the stability diagram of the benzene 
I-SET. 



V. CONCLUSIONS 

In this paper we analyze the transport characteristics 
of a benzene I-SET. Two different setups are considered, 
the para and the meta configuration, depending on the 
position of the leads with respect to the molecule. 

Within an effective p z orbital model, we diagonalize ex- 
actly the Hamiltonian for the molecule. We further apply 
a group theoretical method to classify the many-body 
molecular eigenstates according to their symmetry and 
quasi-angular momentum. With the help of this knowl- 
edge we detect the orbital degeneracy and, in the para 
configuration, we select the states relevant for transport. 

We introduce a generic interference condition (Q} for 
I-SETs in terms of the tunnelling transitions amplitudes 
of degenerate states with respect to the source and drain 
lead. By applying it to the benzene I-SET we predict the 
existence of interference effects in the meta configuration 

In order to study the dynamics of the molecular I- 
SET, we use a density matrix approach which starts from 
the Liouville equation for the total density operator and 
which enables the treatment of quasi-degenerate states. 

The stability diagrams for the two configurations show 
striking differences. In the linear regime a selective con- 
ductance suppression is visible when changing from the 
para to the meta configuration. Only transitions between 
ground states with well defined particle number are af- 
fected by the change in the lead configuration. With the 
help of the group theoretical classification of the states 
we recognize in this effect a fingerprint of the destructive 
interference between orbitally degenerate states. We de- 
rive an analytical formula for the conductance that re- 
produces exactly the numerical result and supports their 
interpretation in terms of interference. Other interfer- 
ence effects are also visible in the non-linear regime where 
they give rise to NDC and current blocking at the border 
of the 6 particle Coulomb diamond as well as to cur- 
rent suppression for transitions between 7 and 6 particle 
states. 

We provide a detailed discussion of the impact of the 
reduced symmetry due to linking groups between the 
molecule and the leads or to an electrostatic potential 
drop over the molecule. We classify different transport 
regimes and set up the limits within which the discussed 
transport features are robust against perturbations. We 
identify in the quasi-degeneracy of the molecular states 
the necessary condition for interference effects. 
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APPENDIX A: GME AND CURRENT IN THE 
NON-SECULAR APPROXIMATION 

The bias and the contact perturbations in our model 
for a benzene I-SET lower the symmetry of the active 
part of the junction and consequently lift the degeneracy 
that appeared so crucial for the interference effects. The 



robustness of the latter relies on the fact that the nec- 
essary condition is rather quasi-degeneracy, expressed by 
the relation SE -C ST. 

Nevertheless, if the perfect degeneracy is violated, the 
secular approximation applied to obtain Eq. Q does not 
capture this softer condition. We report here the general 
expression for the generalized master equation and the 
associated current operator in the Born-Markov approxi- 
mation and under the only further condition (exact in ab- 
sence of superconductors) that coherences between states 
with different particle number are decoupled from the 
populations and vanish exactly in the stationary limit: 



J 



CT EE' 



-(E-E>)a» w + 

E t Pne W 



--p a (F - H&J + f-(F - H° en ) 
L Pa (H° en -F) + f+(Hl n -F) 
Pa (F-H° en )+f-(F-H° en ) 



7T 

i 



r I CT FE' 



-p a (H° cn -F) + f+(H° cn -F) 



4r f FNE' 



arFF' 



E "Tp^NE \ S aT <J? F }d aT 



(Al) 



-p a (E' - F') + f+(E' - F') - - Pa (E -F) + f+(E- F) 



d a N+1 r/ t 

ar FF' "or 



- Pa (F' - E') + f-(F' - E') - - Pa (F - E) + f-(F - E) 

7T 7T 



T'NE' 



where <Tg E / is, differently to Eq. ([7]), in the Schrodinger 
picture. Eq. ((7|) represents a special case of Eq. (|A1[) in 
which all energy spacings between states with the same 
particle number are either zero or much larger than the 
level broadening hT. The problem of a master equation 
in presence of quasi-degenerate states in order to study 
transport through molecules has been recently addressed 
in the work of Schultz et al. 3 ® . The authors claim in 
their work that the singular coupling limit should be used 
in order to derive an equation for the density matrix in 
presence of quasi-degenerate states. Equation (|Aip is 
derived in the weak coupling limit and bridges all the 
regimes as illustrated by Fig. HO 



The current operators associated to the master equa- 
tion just presented read: 



NEFt 



d) 

"cut 



-p a (E-H° cn ) + f-(E-Hl n ) 



-~ Pa (F-H° en ) + f-(F-Hl 



-~ Pa (Hl a -E)+f+(H? 



ben 



E) 



--p a (Hl n -F) + f+(Hl a -F) 



d aT 

(A2) 



where a = L,R indicates the left or right contact. 
Nevertheless, within the limits of derivation of the mas- 
ter equation, this formula can be simplified. Actually, if 
E — F < hT, then F can be safely substituted with E 
in the argument of the principal values and of the Fermi 
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functions, with an error of order ^—S- < -fi^p which is 
negligible (the generalized master equation that we are 
considering is valid for hT <C k^T). The approximation 
E ~ F breaks down only if E~F ~ k^T, but this implies 
E — F 3> ftr which is the regime of validity of the secular 
approximation. Consequently, in this regime, terms with 
E ^ F do not contribute to the average current because 
they vanish in the stationary density matrix. Ultimately 
we can thus reduce the current operators to the simpler 
form: 



la - T a ^ ^Ne| + <4t [fa ( E - H bcn)] d ar 
NEt 

-d aT [f+(H° en -E)] 4 r }, 



(A3) 



which is almost equal to the current operator correspond- 
ing to the secular approximation. The only difference is 
here the absence of the second projector operator that al- 
lows contributions to the current coming from coherences 



between different energy eigenstates. 



APPENDIX B: ANALYTICAL FORMULA FOR 
THE LINEAR CONDUCTANCE INCLUDING 
THE ENERGY NON-CONSERVING TERMS 



In the derivation of the conductance formula (PCJ)1 we 
neglected the energy non-conserving terms in the Eq. ([7]). 
Since in the GME they appear only in the dynamics of the 
coherences between orbitally degenerate states, Eq. ([20)1 
is exact for the para configuration, where the orbital de- 
generacy is cancelled. This is not the case in the meta 
configuration where the orbital (quasi-)degeneracy is es- 
sential for the description of interference. Thus we de- 
rived a generic analytical formula for the conductance, 
taking into account the energy non-conserving terms. It 
reads 



J 



G n ,n+i(A£) = e 2 rA N ,N+i 



SnSn+ifi&E) 



(Sn+i ~ S N )f(AE) + S> 



Here, it is T = Tl = Tr. Ajv,jv+i is the overlap fac- 
tor introduced in Section fill Al Eq. (|21[) . The auxiliary 
function aux(SN, Sn+i) in the correction term is zero if 
there are no orbitally degenerate ground states involved 
in the transition. If one of the participating states is or- 
bitally degenerate it is aux(SN, Sn+i) = 1- The sign 
in f ± (AE) is defined as follows: f + (AE) has to be 
used if the N particle ground state is orbitally degen- 



J 



aux(S , N, S'n+i)3'P 2 



16A 2 , N+1 (S N S N +i) 2 (/ ± (Ai^)) 2 + V 2 



(Bl) 



r 



erate. If instead the N + 1 particle ground state exhibits 
orbital degeneracy, f~(AE) has to be inserted. The 
energy non-conserving terms are included in the factor 
F = T'lJvw^o — ^R-lvbia^o- It is defined only if a de- 
generate state is participating transport. In case that 
e.g. the N particle ground state is orbitally degenerate, 
V a with a = L,R reads 



EE 



E'Z rim 



-p a (E g , N -E') 



EE 

W ,1 nm 



(N-l,E'l\d aT \N e ,n)(N e ,m\di T \N-l,E'l) 



-p a (E'-E K 



(N + 1, E'l\dl T \N g , n) (N g , m\d aT \N + 1, E' I), 



where p a (a;) = -Retp \ + ^{x-n a ) 
digamma function, as defined in Section 



and ip is the 



APPENDIX C: TRANSITION PROBABILITIES 
FOR THE 6 7 TRANSITION 

In the calculation of the overlap factor A in Sec- 
tion UlIA] we used the relation 



\(6 g \d L \7 e ,l = 2)\ 2 = \(6 g \d L \7 g ,l = -2)f 



(CI) 



for the transition probabilities between the 6 particle 
ground state and the 7 particle ground states \7 s ,l), 
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where I is the eigenvalue of the quasi-angular momen- 
tum. This relation is now to be proved. 
Again, we take advantage of the symmetry properties of 
the molecular states with respect to the a v operation and 
to the rotation operator for rotations about a discrete 
angle <j> = ^p, as introduced in Section Ul CI The starting 
point is the generic relation between these two operators: 

7l<p(T v = (JvlZ-c/,. (C2) 

We can now apply both sides of this relation to the 7 
particle ground states \7 g ,l = ±2): 

U^a v \7 g ,l = ±2) = c7 v ft_ |7 g) Z = ±2). (C3) 

The 7 particle ground states \7 g ,l = ±2) are eigenstates 
of each 72.0, and the corresponding eigenvalues are phase 
factors: 

720|7 g , I = ±2} = e^ 2 '^|7 g , I = ±2). (C4) 
Thus, Eq. dCU) becomes 

7^(<7 v |7 g ,Z = ±2)) = e ±24 *(a v \7 s ,l = ±2)). (C5) 



Yet, according to Eq. (|C4|) . this equation can only be 
valid if 

( 7 v |7 g ,Z = ±2) = A|7 g ,Z = T2). (C6) 

and, since cr 2 = 1, A can only be a phase factor. For the 
calculation of the transition probabilities we use further 
the property o~\a v — 1. Since the left contact atom (atom 
1) lays in the reflection plane a v , it is: cr v di,crl = d^. 
Also, since the symmetry of the 6 particle ground state 
is Ai g , it is: <7 v |6 g ) = |6 g ). Under these considerations, 
we can write for the transition probability to the state 
|7 g) Z = 2>: 

|(6 g |d L |7 g ,Z = 2>| 2 = \(6 g \ala v d h ala v \7 s J = 2)\ 2 = 
= \(6 g \d L a v \7 g ,l = 2)\ 2 = 
= |(6 g |d L |7 g ,Z = -2>| 2 . 

(C7) 
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